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This note is about a remarkable law of nature discovered by George Polya [6]. 
Consider a particle situated at a given point of the integer lattice Z d . Suppose that, 
at each tick of the clock, the particle jumps to a randomly selected neighbouring 
lattice point, with equal probability of jumping in any direction. In other words, 
this particle is executing the simple random walk on 1 d . 

A random walk is said to be recurrent if it returns to its initial position with 
probability one. A random walk which is not recurrent is called transient. Polya's 
classic result is the following. 

Theorem 1. The simple random walk on Z d is recurrent in dimensions d = 1,2 
and transient in dimension d > 3. 

This note presents a fairly direct proof of Polya's Theorem which differs from 
both the standard probabilistic proof (see e.g. [H Chapter 4]) and the approach 
via electrical networks [3]. The argument presented here combines basic methods 
from combinatorics (decompositions and generating functions), special functions 
(Bessel function identities), quantum field theory (Borel transform), and asymptotic 
analysis (Laplace's method). Since it illustrates these broadly applicable techniques 
working in tandem to establish a foundational result of probability theory, it is of 
pedagogical value. 

1. Loop decomposition 

Let E denote the event that the simple random walk on Z d returns to its initial 
position, and put p = Prob(-E). For n > 1, let E n be the event that the random 
walk returns to its initial position for the first time after n steps. It is convenient to 
set Eo = 0, corresponding to the fact that the initial position of the random walk 
does not count as a return (if it did, the return probability of any random walk 
would be one). The events E n are mutually exclusive for different values of n, and 

E= □ E n . 

Hence 



P = ^Pn, 

where p n = Prob(i?„). 

A loop on Z d is a walk which begins and ends at a given point. It is convenient 
to consider walks of length zero as loops; such loops are called trivial. A non-trivial 
loop is indecomposable if it is not the concatenation of two non-trivial loops. Choose 
a particular point of Z d , and let £ n denote the number of loops of length n based 
at this point. Let r n denote the number of these which are indecomposable. Note 
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that €o = l while ro = 0. Since any non-trivial loop is the concatenation of an 
indecomposable loop followed by a (possibly trivial) loop, the counts £ n and r n are 
related by 

n 

&n ^ ^ Tk^n — k 
k=0 

for all n > 1. Dividing both sides of this equation by (2d)™, we obtain the relation 

n 

In = ^2 qkPn-k 
k=0 

for all n > 1, where as above p n is the probability that the random walk returns to 
its initial position for the first time after n steps, while q n is the probability that 
the random walk is located at its original position after n steps. 
Introduce the generating functions 

00 00 
P{z) = Y J PnZ n and Q(z) = ^q„z". 

n=0 n=0 

The relation between p n and q n is then equivalent to 

P(z)Q(z) = Q(z) - 1. 

Since p n < q n < 1, each of these series has radius of convergence at least one. Thus 
we may consider P(z),Q(z) as analytic functions defined on the open unit disc in 
the complex plane. The function Q(z) is no n- vanishing for z in the interval [0, 1), 
and hence we have 

Since 

CO 

P(l) = J2Pn=P, 

n=0 

we have 

1 

p=l 



Thus Polya's Theorem is equivalent to the statement that the sum Q(l) diverges 
for d = 1, 2 and converges for d > 3. 

2. Exponential loop generating function 

In order to analyze the sum Q(l), we will obtain a tractable representation of 
the function Q(z). This amounts to finding an expression for the loop generating 
function 

00 

L(z) = J2^nZ n . 
n=0 

Indeed, Q{z) = L(^). 
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While the ordinary generating function L(z) is difficult to analyze directly, the 
exponential loop generating function 

oo n 
n=0 

is quite accessible. This is because any loop on 1 d is a shuffle of loops on Z 1 , and 
products of exponential generating functions correspond to shuffles. This is a basic 
property of exponential generating functions which we will review in the specific 
case at hand. For a general treatment, the reader is referred to [3 Chapter 5]. 

In this paragraph it is important to make the dependence on d explicit, so we 
write tffl for the number of length n loops on Z d and Ed{z) for the exponential 
generating function of this sequence. Let us consider the case d = 2. A loop on Z 2 
is a closed walk which takes unit steps in two directions, horizontal and vertical. 
A length n loop on Z 2 is made up of some number k of horizontal steps together 
with n — k vertical steps. The k horizontal steps constitute a length k loop on Z, 
and the n — k vertical steps constitute a length n — k loop on Z. Thus the number 
of length n loops on Z 2 which take k horizontal and n — k vertical steps is 

y kj k n - k > 

since specifying the times at which the k horizontal steps occur uniquely determines 
the times at which the n — k vertical steps occur. The total number of length n 
loops on Z 2 is therefore 



i — n V / 



r Sj Si- Us" 

This is equivalent to the generating function identity 

E 2 (z) = E 1 (z) 2 . 
The same reasoning applies for any d, and in general we have 

E d (z) = E 1 (z) d . 

3. Bessel functions 
Counting loops in one dimension is easy: 

= /(**)' if ™ = 2fciseven 
1 0, if n is odd 

Indeed, any loop on Z consists of k positive steps and k negative steps for some 
k > 0, and the times at which the positive steps occur determine the times at which 
the negative steps occur. Thus 

~ f 2 k\ z 2k ^ z 2k 

Ei(z) 



^\k)(2k)l £^ kW 



k=0 x ' v ' k=0 

Now a minor miracle occurs: the exponential generating function for lattice walks 
in one dimension is a modified Bessel function of the first kind. 
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The modified Bessel function of the first kind, usually denoted I a (z), is one of 
two linearly independent solutions to the second order differential equation 

d 2 d \ 
z 2 -j^ + z—-{z 2 +a 2 ))F(z)=0, aeC. 

This differential equation is known as the modified Bessel equation; it appears in a 
multitude of physical problems, and was exhaustively studied by nineteenth century 
mathematicians. An excellent reference on this subject is [J Chapter 4]. It is known 
that the modified Bessel function admits both a series representation, 

oo ( z \2k+a 

/«(*) = £■ {l) 



k\T(k 

k=0 y 



and an integral representation, 



I a (z)= e( cose)z (sm6) 2a d9. 
v / TrT(a+ |) Jo 

From the series representation, we see that E\(z) = Io(2z), and hence 

E(z) = I (2z) d . 

4. BOREL TRANSFORM 

We now have a representation of the exponential generating function E(z) of 
loops on Z d in terms of a standard mathematical object, the modified Bessel func- 
tion Iq(z). What we need, however, is a representation of the ordinary loop gener- 
ating function L(z). 

The integral transform 

poo 

(Bf)(z)= / f{tz)e-*dt, 
Jo 

which looks like the Laplace transform of / but with the ^-parameter in the wrong 
place, converts exponential generating functions into ordinary generating functions. 
To see why, write out the Maclaurin series of f(tz), interchange integration and 
summation to obtain 

00 7 n poo 

(B/)( z ) = £/(»)(0)- / ee-'dt, 

and use the fact that 

oo 

re^dt = nl 

This trick is standard issue in quantum field theory, where it is constantly used in 
connection with Borel summation of divergent series, see [SJ §2.3]. 
Applying the Borel transform to our problem, we obtain 

poo 

L(z) = BE(z) = BIo{2z) d = I ^tzfe^dt. 

Jo 

Thus we have arrived at the integral representation 



POLYA'S RANDOM WALK THEOREM 



from which wc obtain 



^)=^ = po(f)V^, 

m= [ io ^ e ~ tdt - 

It now remains to prove that this integral diverges for d = 1,2 and converges for 
d > 3. 

5. The Laplace principle 

The divergence or convergence of Q(l) depends only on the behaviour of the tail 
integral 



N 



M l - ]»-'(//. V >: 0. 



This tail behaviour is in turn determined by the behaviour of the integrand as 
t — > oo. 

In order to estimate the integrand, we invoke the integral representation of the 
modified Bessel function stated above. This yields 



where f(6) = ^cos6>, and we can estimate this integral as t — > oo using a basic 
technique of asymptotic analysis known as Laplace 's method. 

The function f(9) is strictly maximized over the interval [0, tt] at the left endpoint 
9 = 0. Thus the integrand e*-^ 9 ) is exponentially larger at 6 = than at any other 
point of this interval. As t — > oo this effect becomes increasingly exaggerated, so 
much so that the integral "localizes" at 9 = a in the t — > oo limit. To quantify this, 
note that /'(0) = 0, /"(0) < 0, and consider the quadratic Taylor approximation of 

/(*) = 

/(<?)« /(0)-|/"(0)|y. 

Replacing with its quadratic approximation, we obtain the integral approxi- 
mation 

Jo Jo 
Extending the integral on the right over the positive reals results in a half a Gaussian 
integral, which can be computed exactly: 



h °V*l/"(o)l^^ 

'0 

Thus we expect that 



Jo 



2i|/"(0)|' 



o V2*l/"(0)l 
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is an approximation of our integral whose accuracy increases as t — > oo. Laplace's 
principle [2, §5.2] ensures that this is indeed the case: 

Jo \W'(0)|' 
where the notation F(t) ~ G(t), t — > oo means that limt_ ) . 00 = 1. 
Putting everything together, we have the asymptotic formula 

r -t 

IdJ 6 ^ constant • £ 2 , t — > oo. 

Thus we conclude that the recurrence or transience of the simple random walk on 
Z d is equivalent to the divergence or convergence of the integral 

poo 

/ t- d/2 dt, N > 0. 
Jn 

Since this integral diverges for d = 1, 2 and converges for d > 3, Polya's Theorem 
is proved. 
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